Microscopic description of the equality between violation of fluctuation-dissipation 

relation and energy dissipation 
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In systems far from equilibrium, the fluctuation-dissipation relation is violated due to the lack 
of detailed balance. Recently, for a class of Langevin equations, it has been proved that this 
violation is related to energy dissipation as an equality [T. Harada and S. Sasa, Phys. Rev. Lett., 
in press; cond-mat/0502505 . We provide a microscopic description of this equality by studying a 
non-equilibrium colloidal system on the basis of classical mechanics with some physical assumptions. 

PACS numbers: 05.40.-a, 05.20.Gg 
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The construction of non-equilibrium statistical me- 
chanics that is useful for systems far from equilibrium 
is a fundamental problem in theoretical physics. Let us 
recall that linear response theory had been formulated 
for the microscopic description of universal relations es- 
tablished by Einstein, Nyquist, and Onsager pj. Thus, it 
is a significant first step to provide a microscopic descrip- 
tion of formulae derived phenomenologically for systems 
far from equilibrium. 

Recently, an interesting equality relating the violation 
of the fluctuation-dissipation relation (FDR) with energy 
dissipation has been found for Langevin equations under 
non-equilibrium conditions 0, 01- Here, note that the 
FDR is a fundamental relation proved in a linear response 
regime around equilibrium states Q, but it is violated 
in systems far from equilibrium Recent studies have 
revealed that the idea of effective temperature is useful to 
characterize the FDR violation for glassy systems 
and steady state systems 0, 0, 0] . In contrast to these 
studies, this equality claims that the FDR violation is 
characterized by an energetic quantity. 

The equality has been proved for a wide class of 
Langevin systems including many-body systems, time- 
dependent potential systems, and systems in contact with 
many heat reservoirs Q. We then search for such an 
equality in more general non-equilibrium systems that 
are not necessarily described by a Langevin equation. To 
this end, we start with investigating a non-equilibrium 
colloidal system on the basis of classical mechanics. Be- 
cause it is highly expected that the motion of a colloidal 
particle is described by a Langevin equation, the FDR vi- 
olation should be related to the energy dissipation even if 
we describe the system on the basis of classical mechan- 
ics. Here, the energy dissipation is given by the energy 
transfer from the center of mass of the colloidal particle 
to the other mechanical degrees of freedom in the classi- 
cal mechanical description. In this paper, we derive an 
expression of the FDR violation for the classical mechan- 
ical description of the system far from equilibrium. By 
considering the physical conditions of the system, we red- 
erive the equality reported in Ref. Q from the expression 



of the FDR violation. 

Model: Specifically, we study a system of one col- 
loidal particle suspended in a three-dimensional liquid 
confined in a region where —L/2 < y < L/2 and 
—L/2 < z < L/2. Periodic boundary conditions are 
imposed in the x direction. Let Z = (R, P) be the posi- 
tion and momentum of the center of mass of the colloidal 
particle. The colloidal particle is driven by an external 
field E(i) = (Eq + ef p (t), 0,0) and is subject to a spa- 
tially periodic potential Uq(R), where Eq is a constant 
force to realize a non-equilibrium steady state and e/ p (t) 
is a small probe force to investigate the response of the 
system in the steady state. 

The center of mass of the colloidal particle interacts 
with the other mechanical degrees of freedom of the 
colloid-liquid system. Their dynamical degrees of free- 
dom are represented by Y = (ri, pi, • • • , rjy, Pjv)- Fur- 
thermore, we introduce thermostated walls as a net- 
work of boundary particles, whose dynamical degrees of 
freedom are represented by B = (ui, gi, • • • , u/v, gjv')- 
These boundary particles are localized around the bound- 
aries given by y = ±L/2 and z = ±L/2, and interact 
with both the colloidal particle and the molecules in the 
liquid. 

All the interaction potentials including Uq (R) are rep- 
resented by a Hamiltonian H(Z, Y, B). Note that a con- 
stant force Eq cannot be described by the Hamiltonian 
because Eq is regarded as a non-potential force due to 
the periodic boundary condition in the x direction. The 
result presented below does not depend on the details 
of the system. We assume that in the equilibrium case 
(Eq = 0), the statistical properties of a variable set 
(Z, Y, B) are described by a canonical distribution with 
temperature T of the thermostated walls. We also as- 
sume that there exists a steady distribution when Eq ^ 0. 
In order to control the temperature of a finite system, we 
introduce a Nose-Hoover thermostat only for the bound- 
ary particles [ll) . The thermostat may introduce an un- 
physical effect, but we expect that the influence of the 
thermostat vanishes in the limit L,N — > oo with N/L 3 
fixed. Furthermore, when the mathematical rigor is not 
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critically taken into account, the analysis developed be- 
low can be applied to a Hamiltonian system without a 
thermostat in this limit. 

In this mechanical system, the motion of all the parti- 
cles is described by 



where we have defined 
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with P = MdR/dt, p.; = mtdri/dt, and = m^dui/dt, 
where M, mi, and m° are the masses of the correspond- 
ing variables. We express these equations as 



(5) 



where T = (Z,Y,B,X) with A = (Ai, • • • , Ajv)- 

Formal analysis: We consider the time-dependent 
distribution function f(T, t) with an initial condition that 
/(r,to) = /c(T), where 



/c(r) 



1 -/MTo(r)-/3TX£'i A ?/ 2 



(6) 



with i? (r) = ff(r) - f7 (R)- ^ is the normalization 
constant. The distribution function f(T,t) satisfies the 
equation 

df(T t) d d 

= - w ($(r)/(r,t))-e— (*!(t)/(r,t)) . (?) 

Note that / c (r) is the stationary solution of J7J when 
E(i) = and C/ (R) = 0. 
We first set 



/(r,t) = / c (r) e w r ' 4 ). 

Substituting this expression into 0, we obtain 



(8) 



d<KT,t) 
dt 



.* ( r,MM + F( R ).^ 



where F(R) is defined as 

F(R) = E e x 
We solve © formally as 

0(T,t) = / dse A ( s -*) 

«/ to 



ggo(R) 
<9R 



(9) 



(10) 



F(R) 



ago 



+ eJ[s) dP x eJ[s) dP x 



(ii) 



a = $(r) 



<9r" 



(12) 



In the argument below, we fix i and denote the solution 
of JSJ with e = 0, which satisfies T(t) = T, as T(s) with 
to < s < t. That is, T(s) is regarded as a function of I\ 
Then, for an arbitrary function A(T), we obtain 



dA(T(s)) 
ds 



= AA(T')\ r , =r 



Therefore, noting A(T(t)) = A(T), we write 
e A(s -* ) A(r) = A(T(s)), 



(13) 



(14) 



where the right-hand side is regarded as a function of T 
according to the rule mentioned above. Using this ex- 
pression, we rewrite (|llf) as 



<j>(T,t) = J ds [F(R(s)) • V(s) 



dP x (s) 



as) 



dP x (s) 

where we have defined V(s) = P(s)/M. 
Now, let us expand 4>(T, t) as 

^(r, t) = (r, t) + (r, t) + o(e 2 ). (ie) 

Substituting this expression into l|15|) . we arrange the 
terms according to the powers of e. From the terms in- 
dependent of e, we obtain 



with 



<!>W(r,t) = w(r,t ,t;t) 



W(T,t ,s;t) = f dR(s')-F(R(s')), 

■/to 



(17) 



(18) 



where W(T,to,s;t) represents the accumulated work 
done by F during the time interval [ta,s] for the tra- 
jectory satisfying T(t) — T. The expression of the distri- 
bution function © with (|17|) and l|18|) is similar to that 
proposed by Zubarev [l2j and McLennan |l3J. Next, the 
terms proportional to e yield 



^(r,^£, s /p (s) (K( s )-^ r '^^) 



dP x (s) 



(19) 

Furthermore, the average of A(T) by the distribution 
function f(T, t), 



(Mm))*,* 

is expanded in e as 



dTf(T,t)A(T), 



(20) 



(A(T(t))) etto = (A(r(t)))^ + e (A(r(i)))[ o 1) +0(e 2 ). 

(21) 
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Then, from (JSJ, dJ), and {I5J), we obtain 



,£ Ws)(Mt) ™ r »)S ) 



Here, we remark on the steady state of the system. It 
can be expected that (A(T(t))) e t becomes an averaged 
value in the steady state when we take the limit to — * 
—oo. However, f(T,t) itself is divergent in this limit, as 
seen from Ijl7|l and (|18|l . In this paper, we consider the 
limit for the averaged quantities, but we do not study the 
singularity of the distribution function itself. 

Taking this into consideration, setting to ~ * in 
(|22H . we obtain the formula 

t ^-oo\ dP x (s) / to 

(23) 

for t > s. In this formula, the time correlation function 
C(t — s) and the response function R(t — s) are defined 

as 

(V x (t)) ( -l = f d8R(t-8)f>(a), (24) 

J — OO 

C(t-s) = (V x (t)V x (s))^l. (25) 

We also define R(t) = for t < from the causality. 

To this point, no approximation is involved. Then, in 
the equilibrium case (Eq = 0), we derive the FDR 



C(t -s)= TR(t - s) 



(26) 



for t > s. This derivation is straightforward when we 
note that W{T(t),t , s;t) = -U (H(s)) + U (R(to)) in 
this case. Here, we have used 



to-*- oo \ dP x (s) 1 1 



(27) 



to 



On the other hand, when Eq ^ 0, the second term on the 
right-hand side of 123(1 takes a finite (nonzero) value in 
general. Thus, l|23|l provides an expression of the FDR 
violation in the mechanical description. 

Physical consideration: Let us estimate the quantity 
dW(T(t),t ,s;t)/dP x (s) in J23J) for the system consid- 
ered in this study. Our estimation is based on the two 
important assumptions on the time scales. First, let r m 
be the slowest time scale of phenomena that a set of vari- 
ables (y, B, X) exhibits. This time scale r m is expected to 
be much smaller than the relaxation time of the velocity 
of the colloidal particle, which is denoted as Tm. Con- 
sidering this physical expectation, we assume that there 
exists a time scale Ai satisfying 



We choose such a time scale Ai and hereinafter fix it. 
Using this Ai, we define the time-averaged quantity as 

__ ^ i p 2 ds'A(T(s')). (29) 

A l Js-Aj/2 

The second assumption on the time scale of the sys- 
tem is that the relaxation time of the velocity tm is much 
smaller than the typical time scale of the position varia- 
tion r p that is determined by the characteristic length of 
the potential Uo(R). That is, we can choose a time scale 
A2 satisfying 



r M < A 2 < T p . 



(30) 



Now, we apply the perturbation SV x (s) to the phase 
space point T(s). As a result of the perturbation, the tra- 
jectory r(s') with s' < s changes to T(s') + ^r(s'), and 
this change yields the additional work 5W(T(t), to, s; t) 
from qi8[l. It seems reasonable to assume that the re- 
laxation time of the work rate V(s') • F(R(s')) is of the 
order of tm- Thus, from the second assumption H30J1 ■ we 
can estimate 

6W(T(t),t ,s;t)~ f ds'6(V(s') ■ F(R(s'))). (31) 



Here, assuming V • F(r(s')) ~ V(s') • F(R(s')), we write 



V(s') ■ F(R(s')) = V(s') ■ F(R(s') + 77(5'), (32) 

where rj(s) represents the fast part of the work rate fluc- 
tuation whose time scale is much smaller than Ai. Using 
this, we express 



SW(T,t ,s;t)~ / ds , 5[V(s')-F(R( S '))+»7(s')]- 

Js-A 2 

(33) 

From the second assumption (|3()(l . F(R(s')) in the inte- 
grand can be replaced with F(R(s)), and the estimation 



SV(s') = (e ~SV x {s),0) 



(34) 



seems reasonable in the time interval [s — A 2 ,s]. Thus, 
the first term in (|33|l can be rewritten as 

F(R(s)) / ds'8V(Pj ~ F x (R(s))SV x (s)t m . (35) 

Js-Ao 



From this result, we obtain the following expression: 

, dri(s') 



dW(T,t ,s;t) F x {R(s))t m 



dP x (s) 



M 



ds 

a j dP x (s) 



(36) 



T m < Ai < T M - 



(28) 



Recalling the definition of r\ in l|32|l and the first assump- 
tion H28|) . we expect that the second term of (|36|l has 
negligible correlation with V x (t). Further, when we as- 
sume that a friction force for the the colloidal particle is 
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given by —7V, 7 is estimated as M/tm- Using these, the 
substitution of l|36|) into (|23() leads to 



C(t -s) = TR(t -s) + - (V x {t)F x {s)) 
7 



(0) 



(37) 



for t > s. 

Let us define the Fourier transform of C(t — s) as 



dt'e iuit C(t') 



(38) 



Similarly, we define R(u>), while R'(u>) denotes the real 
part of R(uj). Then, from (|37|l . we obtain 



^[C(v)-2TR'(u>)} 

Z7T 



where J = (V^F X )^ is the energy transfer rate from the 
center of mass of the colloidal particle to the other de- 
grees of freedom. This energy transfer rate is interpreted 
as the energy dissipation ratio in the description of the 
center of mass of the colloidal particle. Thus, the equal- 
ity given in l|39[) relates the FDR violation with energy 
dissipation, as presented for a Langevin equation in Ref. 
2J. In this manner, we have rederived the equality on 
the basis of the classical mechanical description. 

Discussion: It might be possible to develop a theory 
for formalizing the above mentioned physical considera- 
tion. In such a theory, the equality given in (GSJ might be 
derived systematically from microscopic dynamics by a 
calculation technique using the separation of time scales. 
The construction of the theory is a future research sub- 
ject. 

Related to this subject, we remark that there exist 
many different expressions of the distribution function. 
For example, when we assume an initial condition in- 
volving H instead of Hq in 10, we obtain a different 
expression of the distribution function. However, in this 
case, we find that it is difficult to relate the obtained 
expression of the FDR violation with the result for the 
Langevin equation. Note that the steady distributions 
for both the initial conditions should be identical when 
the limit to — > —00 is considered. This implies that there 
is an expression that can be treated in a simple manner. 
The clarification of this might provide a key step in the 
systematic derivation of the equality given in (|39|1 . 

A more important question is whether FDR violation 
can be expressed by a form useful for systems in which a 



-J. 



(39) 



Langevin description is not effective. It should be noted 
that one can derive an expression of FDR violation for 
any mechanical system in a similar manner to that used 
for deriving (|23|) . The examples of mechanical systems 
include electric conduction systems, sheared systems, and 
heat conduction systems. However, as mentioned above, 
the obtained expression might have no direct relation 
with measurable energetic quantities. Thus, it is neces- 
sary to find a condition under which the FDR violation 
takes a physically useful form. 

In conclusion, we provide a microscopic description for 
the equality given in (|39|l by analyzing the Hamiltonian 
equation with the Nose-Hoover thermostat at the bound- 
aries. By examining questions arising from this study, we 
wish to obtain a deep understanding of non-equilibrium 
systems. 
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